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Abstract 

The von Neumann entropy density of a block of n spins is proved to be non-zero for 
large n in the non-equilibrium steady state of the XY chain constructed by coupling a 
finite cutout of the chain to the two infinite parts to its left and right which act as thermal 
reservoirs at different temperatures. Moreover, the non-equilibrium density is shown to be 
strictly greater than the density in thermal equilibrium. 
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1 Introduction 

In this letter, we study the large n asymptotic behavior of the von Neumann entropy 

5 W = -tr(p (n) logp (n) ) (1) 

of the reduced density matrix which is the restriction to a subblock of n neighboring spins 
of the non-equilibrium steady state cu + from [7| on the anisotropic XY chain in an external 
magnetic field whose formal Hamiltonian is specified by 

H = ~\ E {C 1 + 1) + " -r) + 2A ^ } • (2) 
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Here, a f , j = 1, 2, 3, denote the Pauli matrices at site x G Z (see © below), the parameter 
7 G ( — 1, 1) describes the anisotropy of the spin-spin coupling, and A G R stands for the 
external magnetic field. In Q (and, for 7 = A = 0, already in H by a different method), the 
non-equilibrium steady state uj + has been constructed in a setting which has become to serve 
as paradigm in non-equilibrium quantum statistical mechanics: a finite cutout of the chain, 
the "small" system, is coupled to the two infinite parts to its left and right acting as thermal 
reservoirs at different temperatures (see SectionEJfor a more detailed description). 
Since the discovery of their ideal thermal conductivity in such states, quasi-one-dimensional 
Heisenberg spin-1/2 systems, made from different materials, have been intensively investigated 
experimentally (see for example K31 the materials SrCu0 2 and Sr 2 Cu0 3 are considered 
to be among the best physical realizations of one-dimensional XYZ Heisenberg models) and 
theoretically (see for example lfT2l l27l l28ll ). Not only this highly unusual transport property 
motivates the theoretical study of correlations in such non-equilibrium models, but the XY 
chain also represents one of the simplest non-trivial testing grounds for the development of 
general ideas in rigorous non-equilibrium quantum statistical mechanics. 

With this motivation in mind, we construct the reduced density matrix of cu + as in ll2"Ull2"6l 
[TBI and prove that the von Neumann entropy <UJ) is asymptotically linear for large block size 
n. Moreover, we show that its non-equilibrium limit density is strictly greater than the limit 
density in thermal equilibrium. 

For further applications of the von Neumann entropy, see Remark 7. 
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Figure 1: The block of n neighboring spins in (UJ). 



2 The non-equilibrium setting for the XY chain 

In this section, we give a brief informal description of our non-equilibrium setting for the XY 
spin model on Z. We refer to J6l [71 [T7| for a precise formulation within the framework of C* 
algebraic quantum statistical mechanics (see also Remark 1 below). 

Consider the XY chain described by the Hamiltonian © and remove the bonds at any two sites. 
Doing so, the initial spin chain divides into a compound of three noninteracting subsystems, a 
left infinite chain TL L , a finite piece Z n , and a right infinite chain 7h R . This configuration is what 
we call the free system whose Hamiltonian 



Hq = Hl + Hn + H 
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is built on the parts Z L , Z n , and Z R according to ©. The infinite pieces Z L and Z R will play the 
role of thermal reservoirs to which the finite system Z n is coupled by means of the perturbation 
H — Hq. In contrast, the configuration described by H, i.e. the original XY chain on the 
whole of Z, is considered to be the perturbed system. In order to construct a non-equilibrium 
steady state cu + in the sense of ll2T1l . we choose the initial state ojq to be composed of thermal 
equilibrium states on the spins on Z^, Z n , and Z R , with inverse temperatures (3 L , f3 D , and f3 R , 
respectively. 

It is well-known that the XY spin model can be described in terms of a model of free fermions 
with annihilation and creation operators b x , b* x , x E Z, satisfying the canonical anticommutation 
relations (CAR), 

{K, by} = 0, {ft*, by} = S X y, (3) 

where {A, B} = AB+BA, the operator b* is the adjoint of b x , and 5 xy is the Kronecker symbol. 
This description is achieved with the help of the Jordan- Wigner transformation |[T9l . 

b x = T^\a{ x) -ia^)/2, (4) 

where ir^ = . . . <J^~^ for x > 1, = 1, and tt^ = . . . for x < 1 (the spin 

(x) 

(jj at site x, j = 1, 2, 3, is given by ... ® 1 2 <S> I2 <8> Oj ® I2 <8> l2--- 5 and 1 = ... <g> 1 2 ® I2 ® 
l m denotes the identity matrix on C m ). Here, a , a 3 is the Pauli basis of C 2x2 , 
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(5) 



Moreover, T has the properties T* = T,T 2 = 1, and Taf'T = 0-{af>), where, for j = 1, 2, 
is defined by 6_{af ] ) = -af ] if x < and 9_{af } ) = } if x > 1. 

Remark 1 The element T stems from the C* crossed product extension by Z 2 of the C* algebra 
& of the spins on Z for the two-sided chain (see |7|; for the one-sided chain, T can be 
dropped): The C* algebra O generated by & and the element T contains the CAR algebra 
A generated by the fermions b x , b* x on Z as a C* subalgebra. The extension to O of the *- 
automorphisms 6 on 6, defined as the rotation of the spins about the z-axis with angle n, 
9{aY > ) = ~' yt f S> f° r j = 1) 2 and^fa^) = a^ x \ yields the decomposition of & and A into even 
and odd parts with respect to parity, 6(A) = ±A, A e O. The even parts coincide, 6+ = A+, 
and form a C* algebra (for the odd parts one has (5_ = TA- and (5_ is merely a Banach 
subspace). Since the extension to O of the "-automorphism groups r and r generated by H 
and H leave the even parts invariant, and since the unique KMS states on & and A coincide 
on & + = A + , the even parts only are of importance in the construction of the non-equilibrium 
steady state u + in 1171. 

Using Kato-Birman scattering theory for the construction of the wave operators on the 1 -particle 
Hilbert space of the Jordan- Wigner fermions Q, the non-equilibrium steady state uj + with 
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respect to the initial state lu and the time evolution r* generated by H has been constructed in 
as the limit 

LU+(b x ) = lim u} (T\b x )) =u o j + (b x ), 
t— > +00 

where 7+ denotes the algebraic analogon of the wave operator on Hilbert space, the so-called 
M0ller morphism on A. Moreover, it has been shown in J71 that cu + is a quasi-free state charac- 
terized by its 2-point operator S (see Appendix lAll. 

u + (B*(f)B(g)) = (f,Sg), (6) 

where, with f) = l 2 (Z), f, gGh©h~h®C 2 (and (/, g) denotes the scalar product in f) © h). 
Here, B is the linear mapping 

/ = (/+, /_) » B(f) = (/+(*) K + /-(x) b x ) (7) 

introduced in Q in the framework of self-dual CAR algebras (see Appendix lAl. It has the 
properties 

{£*(/), B{g)} = (/, g), B*(f) = B(Jf), (8) 

where J is the antiunitary involution defined by (/+, /_) t— > (/_, /+) (the bar denotes complex 
conjugation). Finally, S has been explicitly computed in 0. In the Fourier picture, Z 2 (Z)®C 2 ~ 
L 2 (T) ® C 2 (with T = {z e C I \z\ = 1}), 5 reads 

s (0 = (l + e^ ft «) +ifc «))" 1 J 
with the temperature parameters /3 and 5 given by 

= i(fo + &), <*=~(fti-&)- (9) 
The 1 -particle operators h and have the form 

K€) = (cos£- A) <g> 0-3 - 7Sin£ <g> a 2 , fc(£) = sign(/«(£)) <8> <7 , 
where the functions and are defined by 

«(£) = 2A sin e - (1 - 7 2 ) sin 2£, MO = ((cos £ - A) 2 + 7 2 sin 2 £) 1/2 . 

Expanding s(£) with respect to a , cr 3 from ©, we find the 0-th component of s(£) = s (C)® 
°"o + Z)fe=i ® to look like 

so(0 = 2 ~ ^^(O signK(0, (10) 
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whereas Sfc(£)> k = 1, 2, 3, has the form 

= "^^(0^(0, r(0 = (0,-7 sine, cose -A). (11) 
Here, we used the definition 

= ^(StSUk)) - a,a ' eR - <12) 

In order to study the von Neumann entropy we restrict cj + to spins on the finite configuration 
subset A„ = {l,...,n}, n G N, of the entire spin chain on Z. Due to the Jordan-Wigner 
transformation @, this is equivalent to restricting uj + to the fermions 6 f , 6* at the same sites A n 
(on the even part of the algebra, see Remark 1). Let h n = / 2 (A n ) ~ C n be the state space over 
A n (where here and in the following, h n is considered as a subspace of f) by the trivial injection). 
Since the Fock space #(f) n ) over the n-dimensional state space h n is 2 n -dimensional, the CAR 
algebra A n = A{\) n ) over h n is 2 2n -dimensional, *4 n ~ £(jF(h n )) ~ c2 n x2 n ^ see f or exani pj e 
|[6l or ifTDl p. 15]; we denote by C(H) the bounded operators on the Hilbert space H). Therefore, 
the restriction of u> + to the spins on A n is described by a density matrix G £(jF(h n )) (see 
for example JlOl p. 267]), 

w + (A)=tT(j>WA), AeA n . (13) 



3 The asymptotics of the von Neumann entropy 

In Theorem our main result, we prove that the large n asymptotics of the von Neumann 
entropy from ^ in the non-equilibrium steady state uj + is, to leading order in n, pro- 
portional to the volume of the spin block, and we derive an explicit expression for the non- 
vanishing proportionality constant. To do so, we express with the help of the eigenvalues 
of a Toeplitz matrix which we will construct next. Let us first define the Majorana correlation 
matrix Q n G C 2nx2n by 

(Q n ) kl = u + (B(wW)B(w®)), k, I = 1, 2n, (14) 

where B{f) is given in ©, and = {w+\ io_ ) G h„ © h n is specified, for j = 1, n, by 

w W(j) = (©» fci2i _ 1 = 5 fe , 2 ,_ 1 + ^ fe , 2j . 

Moreover, the matrix r G C 2x2 reads 
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Remark 2 By virtue of ©, the 2n operators d k = B(w^) are Majorana operators, i.e. 

d* k = 4, {4, di} = 25 kl . (16) 



In the following, we define a suitable block Toeplitz matrix T n [a\. The reader not familiar with 
Toeplitz theory may consult Appendix |B] where the notation and the relevant facts are given. 
Moreover, in the proofs below, we will refer to Appendix |A| for an elementary construction of a 
suitable fermionic basis. 

Lemma 1 The imaginary part of the Majorana correlation matrix Q n from (1141) is a skew- 
symmetric 2x2 block Toeplitz matrix T n [a] with symbol a £ L^ 2 , 



2 

tt n = l 2n + iT n [a], a(£) = - 

i 



s 2 (0 + i« 3 (0 s (0-| 



(17) 



where the functions s , s 3 are defined in (fl~Qb and (fTTT) . 



Proof With the Fourier transform F : Z 2 (Z) -> L 2 (T) in the form = £) x6Z /(x) e ia * we 

have Fw^' 1 ' = e k , Fw± k ^ = =Re fc , and S = F*®1 2 s F®1 2 , where we define e x (£) = e lx ^. 
Hence, with 7)± = (1, ±1) £ C 2 , we compute 

3 

= 2 (ej, (s + si) e k ) 

= S jk - f*7^<Pw{£) e-*^. (18) 

In the last equality only, we used the explicit expressions of s and si from (fTOl) and (fTTT) (the 
function is defined in ()12|> ). In the same way, we obtain 



(fi n ) 2 j-i,2fc = 2 (ej, (s 2 - is 3 ) e fc 

2tt //(£) 

(0)2^-1 = 2(e jl (s 2 + is 3 )e k 



m(0 e- i(fc - J)f , 



where we define 



(^n)2j,2fc = 2(e i ,(s -s 1 )e fc ), (19) 



g(0 = -7sin{ - z(cos{ - A). 
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Hence, the matrix T n [a] = (Q n — is a block Toeplitz matrix with 2x2 blocks, i.e. 

{£l n )k+2i,i+2i = (Q n )ki- Its symbol a E L?? x2 can be read off from (H5l and (THJt. 

a ^ = i wM) sign«(0 J§ 

Moreover, using the properties (fToT) of the Majorana operators for i(T n [a]) k i = uj + (d k d{) — 8 kh 
we see that T n [a] is real and skew-symmetric, i.e. (A* denotes the transpose of the matrix A) 

T n [a}eR 2nx2n , Tniaf = -T n [a}. 

This is the claim of Lemma d □ 



In the next lemma, we choose a special set of fermions in A n stemming from the block di- 
agonalization of T n [a] by means of real orthogonal transformation V E 0(2n) = {V E 

R 2nx2n | yyt = 

Lemma 2 There exists a set of fermions Ci = of 1 ' E A n , and numbers \\ E R, i = 1, n, 
such that the reduced density matrix p- n ' from (1131) has the form 

P (n) = f[ f^ 1 ^ c* Ci + . (20) 

Proof Let q, i = 1, n, be arbitrary fermions in *4 n and define operators e^L i = 1, n, 

a, (3 = 1, 2, by 

(i) * (j) * (j) (i) * cyw 

e ll ~~ C i C i; e l2 ~~ C i ' e 21 ~~ C «' e 22 ~~ C * C j ■ V^-U 

Since Lemma|51in AppendixlAl states that the 2 2n operators n™=i e a;ft' a i' A.) •••> a m Ai = 1> 2, 
constitute an orthonormal basis in £(#((}„)), we can write p( n ) G £(-5 r (f) n )) as 



P 



{ " ] - e «♦ its* ikv < 2 » 



«i,/3i,...,a„,/3 n =l,2 \ \i=l / / j=l 



The goal of the following is to make a special choice for the fermions c, such that the density 
matrix from (l22l can be written in the form d20l) . For this purpose, let V E 0(2n) be any 
real orthogonal 2n x 2n matrix and choose 



Ci = B(v {i) ), 1 = 1,. ..,71, 



(23) 
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where uW = (v® , v_) G h n © h n is specified, for j = 1, n, by 

v%) = m^V^r))^.^ (24) 

and r is defined in (1131) (note that for the choice V = l 2n we recover the Jordan-Wigner fermions 
hi from ©). Using 

2n 

Ci = E((©^" 1 )^)»-i lfc *, (25) 
fe=i 

the Majorana relations (IToT) . the fact that V G 0(2n), and the properties of r, we see that the 
operators c, are fermionic, 



{q, c,-} = 2 ((©^- 1 )yy*(©^r- 1 *)) 2i _ 1 , 2i _ 1 = (©>i) 2j -i, 2i -i = 0, 
{c*, Cj } = 2(((B n 1 T- 1 )VV t (® n 1 T- 1 %^ 1 ^ 1 = (e>o) 2j --i,2i-i = 5i 



(note that, in general, the CAR are not satisfied if we choose V to be complex unitary or complex 
orthogonal). Next, let us make a special choice for V. Since we know from Lemmadthat T n [a] 
is real and skew-symmetric, there exists a matrix Q G 0(2n) which transforms T n [a\ into its 
real canonical form (see for example [16, p. 107]), 

QT n [a) Q* = ®] =1 (\f ia 2 ) = diag A^) ® za 2 , (26) 

where zHAj"% A^ G R, are the eigenvalues of T„[a]. Hence, setting = Q in (l24h . and using 
d25t and d26l) . we find, similarly as for the CAR above, 

uj + (ciCj) = -(e^xiax- \ k ia 2 )) 2i _ lj2j _ 1 = 0, (27) 
1 1 + A^ 

Now, let us write p( n ) from d22l with the help of this special choice of fermions. Since d27t and 
hold, we have, from Lemma|5]in Appendix|Al that 



iKk Hn^^(&- (29) 



,1=1 / / 8=1 



Therefore, plugging d2^1) into (l2"2l and using (OTT) . dTTl . and (l2"%b . we find that the reduced 
density matrix p^ n > takes the form 

n n / 1 + A (n) 1 - A (n) 

p (n) = n (^+( c * c *) + ^+( c i c *) c ^*) = n — ^~ ^ + — ^~ qi 

i=l i=l V 
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which is CD). D 

Next, we formulate our main assertion about the asymptotic behavior of the von Neumann 
entropy <S^ = — tr (p^ log p^) in the non-equilibrium steady state u + characterized by © 
with its reduced density matrix p^ from (fT3l) . The temperatures (3 and 5 are defined in ©. 



Theorem 3 Let uj + be the non-equilibrium steady state of the XY chain with reduced density 
matrix p^ n \ Then, for the temperatures < 5 < (3 < oo, the anisotropy 7 G (—1,1), and the 
magnetic field AgI, the von Neumann entropy S^- 71 ' is asymptotically linear for n — > 00, 

=Cn + o(n) with C = C PAl ,\ > 0, 

and C is given in (l34l) . 

Remark 3 The existence of the limit lim n ^ 00 S^/n for translation invariant states of spin 
systems follows from the strong subadditivity property of the mapping [1, n] — > S^ n \ see for 
example [ 10, p. 287]. In (l34l) . we derive an explicit expression for this limit and show that it is 
strictly positive. 

Proof The properties d47t and (l48l) from Appendix |B] imply that the spectral radius of the 
Toeplitz matrix T n [a] is bounded by 

\\T n [a}\\ < esssup ||a(0l|£(c 2 ) = ess sup + fs,^)) = Q, (30) 

f e [0,2tt] £e [0,2tt] 

where, due to the finiteness of the temperatures, g is strictly smaller than 1, 

g = th(/3 fl (l + ]A])/2) < 1. (31) 
The form d20l immediately yields the spectral representation of p^ 1 ' with the 2™ eigenvalues 

- 1 4. C_iV<\( n ) 

A ei ,..., en = n 2 ^ > e 1; ...,e n G{0,l}. (32) 

i=l 

Due to <|30|) and (|3T|> . we have < A eii ..., e „ < 1. Using (l32|) . trp (rt) = 1, and the entropy 
function h : (—1, 1) — > IR + (for which = r/((l + x)/2) with r)(x) = — xlogx — (1 — 
x) log(l — 2) the Shannon entropy; see Figure 2 for h(x)), 

ftW = -i±£lo g fi±^V 1 ^'ogf i ^y (33) 
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Figure 2: The entropy function h(x) from (CHI . 



the von Neumann entropy = —ti(p^ log p^) can be written in the form 

n 

S {n) = ~ E ^,...,e„logA ei ,..., e „ = EMAS n) ). 

ei,...,e„e{0,l} i=l 

To determine the asymptotic behavior of S^ n \ we want to make use of Szego's first limit theo- 
rem for the case of block Toeplitz matrices (see @ p. 202] and d49t in Appendix|B]). To do so, 
we note that the symbol ia G L?? x2 is self-adjoint (with a from (fTTl) ) and that ±Aj G (— £>, £?) 
are the eigenvalues of the self-adjoint Toeplitz matrix zT n [a] = T n [ia] (see (l2"6T) ). Moreover, 
with the definition h(x) = h(x) for x G (—1, 1) and h(x) = for x G R \ (—1, 1), we have 
/i G Co(M) (the continuous functions with compact support) and hence, due to h(—x) = h(x), 
Szego's first limit theorem d49l> implies 



n— >oo ?T, 

f27T 



= 5is tr * (wK)) 
i r 2n df 

= 2 /o ^^ 5 ^ + ^(0) + M^(0-^,/3(0)) (34) 

> > o. 

This proves our theorem. □ 



Remark 4 Contrary to what one would expect naively, the singular nature of the symbol does 
not affect the leading order of the asymptotics of the entropy density: As soon as there is a 
strictly positive temperature in the system, the spectral radius of the Toeplitz operator contracts 
to a value strictly less than 1. The same phenomenon has been observed for the transversal 
spin-spin correlations in this model, see 0. 

Remark 5 In the case 5 = 0, the symbol of the Toeplitz operator becomes smooth. Then, 
second order trace formulas imply that the subleading term has the form o(n) = const + o(l), 
seeElp.133,207]. 
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Remark 6 Note that for (3 — > oo, the right hand side of (l34b vanishes in agreement with the 
findings in the literature. In this case, the logarithmic behavior of the entropy density has been 
derived using the asymptotics of Toeplitz operators with Fisher-Hartwig symbols, see lfT8l . 

Remark 7 The von Neumann entropy is widely being used as a measure of entanglement 
in the ground state of a variety of quantum mecanical systems. In ll20l l26l fl~8l for example, it 
has been shown that the entanglement entropy tends to a finite saturation value, whereas at a 
quantum phase transition (see 11221 s ). it grows logarithmically with the size n of the subsystem. 
These results can be related to conformal field theory computations of the geometric entropy, 
see, for example, El EH M O- 

Corollary 4 For large n, the non-equilibrium entropy density is strictly greater than the entropy 
density in thermal equilibrium, i.e. Cp^x > Cp,o,y,x- 

Proof From (l34T> we have 



Remark 8 In lfT3ll . following [1 ], a non-equilibrium steady state has been constructed by im- 
posing a current on the system (these states differ in general from the non-equilibrium steady 
states in the sense of lETL see also 0). The content of the corollary is in agreement with [ 13] 
in the sense that the von Neumann entropy increases in the current carrying state (by doubling 
the prefactor of the asymptotically logarithmic behavior). 

A Quasi-free states on self-dual CAR algebras 

Quasi-free states on self-dual CAR algebras have been introduced in Q. We briefly review 
these notions. A self-dual CAR algebra over a Hilbert space K, with scalar product (•, •) and 
antiunitary involution J is a C* completed *-algebra generated by B(f), f E /C, B*(f), and an 
identity 1 which satisfy 




(/i(th(/3 R/ i(0/2)) + fc(th(/?flM0/2 - Sp(0))) > CW- 



□ 



{B*(f),B(g)} = (f,g) t B*(f) = B(Jf) 



(35) 




(36) 



m 
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(B(f^)...B(f^)) = ^sign7rn^(^(/ (7r(2i_1)) )^(/ (7r(2i)) )), 



(37) 



TT i = l 
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for m G N, and the sum runs over all permutations n in the permutation group 5 2m with 
signature sign7r which satisfy ix(2i — 1) < 7r(2«) , 7r(2i + 1) (the pairings of {1, 2m}). Such 
a state is characterized through 

u(B*(f)B(g)) = (f,Sg) 
by its 2-point operator S on /C with the properties (see 0) 

0<5 , <1, S + JSJ=1. 

Remark 9 The definition of B(f) in © and of J after ®, with / G /C = h © h, is a special 
case of a self-dual CAR algebra, and the non-equilibrium steady state lo + from © is a quasifree 
state in the sense (EH), (OH), see 0. 

Let q, z = 1, n, be any fermions in A n ~ ~ C 2 ™ x2 '\ 

{q, cj = 0, {c*, cj = 4,-, (38) 



and define operators e^L % — 1, n, a, f3 — 1, 2, by 



(0 * « * (0 00 

-11 — e 12 ~~ c i> e 21 ~~ c «5 e 22 



(39) 



It follows from the CAR d38t that the operators e^L have the properties 



= (-l) (Q+/3)(7+5) eS } eS, for i * j, (40) 
= <W2- (41) 



Remark 10 Note that, due to (f35b. the operators q = B(f^) are fermions d38t . if and only 
if (Jf {i \f ij) ) = and /&')) = In the case at hand, the fermions are given in (IziD 
(there is a ^isomorphism between the CAR algebra A over 1} and the self-dual CAR algebra 
over fj © h, see Q). 

Lemma 5 Le? e^l be defined by (l39l) for fermions Ci = B(f^), and let on be a quasi-free state 
on A n which satisfies 

Lo(ciCj) = 0, uj(c*Cj) = 5ijUj(c*Ci), i,j = 1, ...,n. (42) 
Then, to factorizes on fllLi e 2ft as 



to 



,i=l / i=l 
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Proof Assume first that the number of £>(■) in YYi=i e a a a * s °dd. Then, there exists an i £ 
{1, n} such that oti ^ Pi and, hence, due to d36l) . (l43l) holds. Now, let the number of £>(•) in 
riLi e a- ( 3 i be even. If oii 7^ $ for some i G {1, n}, then every term in the sum (l37t contains 
a factor of the form u(B(J ei f^)B(J € ^f^)) with z 7^ j and some e i: ej G {0, 1}. Therefore, 
due to (l42l . we have 

(n \ n / n \ 

i=i / i=i \j=i / 

For the same reason, applying (l37t to the right hand side of d44b . the only nonvanishing term is 
the one for which n is the identity permutation. This implies (l43l . □ 



nn-rnvrm 

Figure 3: The pairings of the right hand side of d44l) for n = 2. 

Lemma 6 Le? e^l £e defined by d39t /"or fermions Ci = B(f^). Then, the set of operators 

lir=i e «ift' "bAi-i^n,^ = 1,2, w an orthonormal basis of the Hilbert space An with 
scalar product (A, B) 1— > tr(A*B). 

Proof We show that the 2 2n vectors Yl7=i e a-ft w ^ m a i> A' •••> a «' ^ = 1; 2, are orthonormal 
with respect to the scalar product *4. n x *4.„ 9 (A, 5) i-> tr(A*i?) by proving 

n^iKU =nww (45) 

i=i \j=i / / i=i 

For this purpose, we push the j-th term in the second product on the left hand side of d45l) 
to the left using (l40b until it it hits the j-th term in the first product and apply (BIT) (note that 

(7)* (i) 

I s - = e £ yP" Doing so ' from i = n un til 3 — 1> we arrive at 

71 / 71 \ ft 71 — 1 

nc, iK4 = n**» nM> ( " +i ' ,E ^ , "*' +v irv <«> 

i=l \i=l / i=l 3=1 Z=l 

Since tr(-)/2 n is a quasi-free state on *4. n (see for example Q) which, due to the CAR d38t and 
the cyclicity of the trace, satisfyies (l42l) . we have, using (l43l) and (l46b . 

(n / n \ *\ n n— 1 n 

j=i \j=i / / i=i j=i ;=i 
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(of course, (l43t can easily be verified directly for tr(-)/2 n ). Finally, using again the CAR d38l) 



and the cyclicity of the trace, we have tr(e^ Qi )/2 n = 1/2 and d45t follows 



□ 



B Toeplitz operators 

Toeplitz operators p. 185] Let IVgN. We define the space l 2 N of all C^-valued sequences 

/ = {/i}r=i> fi G C", by 



1/2 



& = {/ : N 



< oo}, 



En/- 



i\\C N 



1=1 



where || ■ || c jv denotes the Euclidean norm on C^. For N = 1 we write I 2 = l\. 

Let {a x } xe z be a sequence of N x N matrices, a x £ C NxN . The Toeplitz operator defined 
through its action on elements of l 2 N by / i— > {J^ili a «-j * s a bounded operator on l 2 N , if 

and only if 



o 



2tt 



for some a 6 L^ xN (see @ p. 186]), where we define (with T = {z £ C | |z| = 1}) 

L™ xN = {0 : T - C* xJV | ^ £ L~(T), 2, J = 1, N}. 
In this case, we write the Toeplitz operator as 



T[a] 



a a_i a_ 2 ... 

(l2 dl Oq 



The function a £ L^ xN is called the symbol of T[a\. If N = 1 the symbol a £ L°° = Lf xl 
and the Toeplitz operator T[a] are called scalar, whereas for N > 1 they are called block. 

For n £ N, let P n be the projections on l 2 N , 

Pn{{fl, •••) /n, /n+1, •••}) = {/l, /n> 0, 0, ...}. 

With the help of these P n , we define the truncated Nn x iVn Toeplitz matrices as 

T n [a] = P n T[a]P n \ taaPa , 
where ran A denotes the range of the operator A. 
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Norm of Toeplitz operators [ 9 , p. 1 86] The norm of a Toeplitz operator is related to its symbol, 

Mil = Halloo, (47) 

where || a|| ^ is defined to be the operator norm of the multiplication operator acting on ©^L 2 (T) 
by multiplication with the matrix function a £ L^ xN . We have 

Halloo = esssup IKOIU^)' ( 48 ) 

€G[0,2vr] 

where || ■ H/vow) is the operator norm induced by the I 2 norm on C N , see also [8, p. 101]. 

Szego's first limit theorem in the block case p. 202] Let a £ be self-adjoint, a* = 

a, and let \[ , be the eigenvalues of the block Toeplitz matrix T n [a\. Then, for any 
continuous function / with compact support, / £ C (M), we have 

lim TtEM" 1 ) = \ ? tr /K0)- (49) 

n.->oo JSn I\ J 2tt 
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